Fluctuations of company yearly profits versus scaled revenue: 
Fat tail distribution of Levy type 

H. EDUARDO ROMAN^ , RiCCARDO A. SlLIPRANDI^ , CHRISTIAN DOSE^ , CLAUDIA RlCCARDI^ AND 

Markus Porto^ 

^ Dipartimento di Fisica, Universitd di Milano-Bicocca, Piazza delta Scienza 3, 20126 Milano, Italy 
^ Hewlett-Packard, Via Giuseppe Di Vittorio 9, 20063 Cernusco sul Naviglio (MI), Italy 

^ Institut fiir Festkorperphysik, Technische Universitat Darmstadt, Hochschulstr. 8, 64289 Darmstadt, Germany 



PACS 89.65.Gh - Economics; econophysics, financial markets, business and management 
PACS 05 . 45 . Tp - Time series analysis 

PACS 05 . 40 . -a - Fluctuation phenomena, random processes, noise, and Brownian motion 

Abstract. - We analyze annual revenues and earnings data for the 500 largest-revenue U.S. 
companies during the period 1954-2007. We find that mean year profits are proportional to 
mean year revenues, exception made for few anomalous years, from which we postulate a linear 
relation between company expected mean profit and revenue. Mean annual revenues are used to 
scale both company profits and revenues. Annual profit fluctuations are obtained as difference 
between actual annual profit and its expected mean value, scaled by a power of the revenue to 
get a stationary behavior as a function of revenue. We find that profit fluctuations are broadly 
distributed having approximate power-law tails with a Levy-type exponent a ~ 1.7, from which 
we derive the associated break-even probability distribution. The predictions are compared with 
empirical data. 



Predicting forthcoming year company profit is difficult 
due to the many unknown variables determining the actual 
earnings scenario. This intrinsic uncertainty in economy's 
evolution makes earnings forecasts not to be correlated 
to actual earnings with the desired accuracy. A conse- 
quence is that, often, stocks with highest earnings fore- 
casts dramatically underperfom those with poor forecasts 
(see e.g. [1]). 

Indeed, company earnings may undergo dramatic vari- 
ations from year-to-year, even over shorter time scales, 
leading to huge movements in public company stock (see 
e.g. [2]). A less volatile quantity is represented by to- 
tal company revenue, but also in this case revenue vari- 
ations may yield conspicuous changes in the underlying 
stock price. Interestingly, the connection between stock 
price (i.e. market value) and revenue is still surrounded 
by many open questions which are awating for further re- 
search (see e.g. [3]). Clearly, the question arises of how to 
estimate in a more realistic fashion profit fluctuations, and 
therefore attempting to improve the accuracy of earnings 
predictions, the latter being closely related to the issue 
of profitability or break-even point [4]. Several attempts 
have been made in order to incorporate a stochastic be- 



havior of profits into the analysis (see e.g. [5-10]), in which 
fluctuations are assumed to be normally distributed. 

From a fundamental point of view, one may wonder 
whether the above difficulties can be mitigated to some 
extent by modifying the way the problem is approached. 
In physical many-body systems for example, a flrst, real- 
istic solution to a problem can be achieved if one resorts 
to the so-called mean-field approximation, in which a sin- 
gle particle 'sees' an average field due to the remaining 
particles in the system (see e.g. [11,12]). Particle-particle 
correlations and fluctuations of physical quantities can be 
incorporated into the formalism at a later stage once the 
mean-field solution of the problem is known (see e.g. [13]). 
How can we apply this idea to the study of profit fluctu- 
ations of real companies, which can be viewed as a many- 
particle system of interacting economic units? Is it pos- 
sible to come up with the strong fluctuations observed in 
company profits? 

Profit fiuctuations can be naively evaluated by looking 
at their relative variations say, from year to year. As a 
matter of fact, however, profit is closely related to revenue 
and production costs (see below) and a different approach 
based on these interrelations could be explored. 
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In this Letter, we suggest that revenue can be taken 
as the independent, driving variable and present a novel 
analysis of profit fluctuations based on this assumption. 
We support our premises using market data from U.S. 
companies on an year-to-year basis over a period of 54 
years. The analysis of the empirical data suggests a form 
for the expected mean profit, being a function of com- 
pany revenue, with respect to which earnings fluctuations 
can be determined. The latter turn out to be dependent 
on revenue, suggesting that they are not stationary as a 
function of revenue. Invoking then the condition of sta- 
tionarity, the fluctuations are scaled by a power of revenue 
with an exponent 77 in the range 1/2 < 77 < 1. The proba- 
bility distribution function of scaled fluctuations displays 
slowly-decaying tails which turn out to be of Levy type. 
A further analysis on the resulting break-even point yields 
a prediction, supported by empirical data, for the prob- 
ability of profitability, enlightening the role that market 
fluctuations play in the problem. 

In the following, we briefly review the standard cost- 
volume-profit (CVP) analysis, from which we derive our 
main conjectures regarding profit fluctuations. Let us con- 
sider a generic (typical) company. According to standard 
CVP analysis [14, 15], we write the profit P as the differ- 
ence between total revenue R and costs, the latter being 
the sum of variable costs 14 and additional (sometimes 
referred to as fixed) costs F > 0, that is 

P = R-{V, + F). (1) 

Further, we write total revenue as R — VsUg, where Vg 
is the sale price of sold unit and rig the total number of 
sold units. Similarly, total variable costs are written as 
Vc = VcJic, where Vc is the cost of produced unit and Uc is 
the total number of produced units. 

In what follows, we assume linear relations between sold 
unit values and produced ones, according to 

Vs ~ as Vc, 1 < Qfs, 

n, ~ 13, nc, < A < 1. (2) 

The above assumed linearities do not preclude the coeffi- 
cients (as,/3s) to be time-dependent, similarly to a sort 
of piecewise linear approximation applied in non-linear 
CVP analysis [16]. Using these relations in Eq. (1), with 
K = R/{asPs), yields 

P = 7.i?-F, with 7s = ^i^^. (3) 

Since we are interested in the typical behavior of compa- 
nies, we write down the above P-R-F relation in terms of 
its mean values, Pq = (P), Rq = (R) and Fq = (F), rep- 
resenting averages of P, R and F over several companies 
for a given time horizon, say a year, i.e. 

Po = {ls)Ro - Fo. (4) 
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Fig. 1: (color online) (a) Mean yearly revenue, Rq [Billions] 
(circles) and mean yearly profit, Po [Billions] (squares) of 
the 500 largest revenue U.S. companies [17] as a function 
of the year. The thick straight line (top) is a fit with the 
form: Ro = ^oexp[(year - 1954)/Po] B, with Ao = 0.27 
and Bo ~ 12. The dashed line (bottom) is the form: Po = 
yliexp[(year - 1954)/Pi] B, with Ai = 0.014 and Bi = 12. 
The thin straight line is an exponential regression for Po over 
the whole years, yielding: Pq = ^2 exp[(year — 1954)/P2] B, 
with A2 — 0.016 and B2 = 13.8, which is not proportional 
to Po- (b) Mean yearly profit, Po, vs mean yearly revenue, 
Po (circles) (from (a)). The straight line is the linear form 
Po = O.O52P0. Anomalous-years results are indicated by the 
down triangles. 

To test this relation, we consider the set of 500 largest 
revenue companies in the U.S. during the period (1954- 
2007) [17], for which yearly values of P and R are avail- 
able. For each year in the database, we calculate the mean 
values Po and Rq in billions (B) of U.S. dollars. Empirical 
results for Rq and Pq are plotted in Fig. 1(a) as a function 
of year. These results suggest that both quantities grow 
exponentially, and that Pq oc Rq if the anomalous years 
1991/92/93 and 2001/02 are excluded. The exponential 
dependence of mean profit and revenue also reflects the 
growth of companies [18], displaying other interesting fea- 
tures described by exponential distribution functions. 

Now, to better appreciate the apparent proportionality 
between mean profits and revenues, we have plotted them 
in Fig. 1(b), suggesting that indeed 

^^0 ^ {l,)Ro, (5) 

with (7g) = 0.052. Slight deviations from linearity can be 
observed in Fig. 1(b) at large revenues, Rq > 10 B, corre- 
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spending to recent last years. This is an indication that, 
possibly, non-linear corrections to the result Eq. (5) are 
playing a role. In keeping with our mean-field strategy, 
however, we will consider such deviations as due to typi- 
cal market fluctuations. Within this scenario, the model 
seems to be consistent with downward profit corrections 
for 2008, and possibly for the next few years, respond- 
ing to a sort of reverse to the mean linear behavior ob- 
tained in Fig. 1(b). As a matter of fact, there is already 
a widespread consensus that 2008 is manifesting a rather 
weak economic environment. Furthermore, the linear re- 
lation between Po and then suggests that 



Fo ~ (7c)i?o, 



(6) 



such that 

(7g)-(7s)-(7c). (7) 

Although we do not know ^7c) explicitely, one can argue 
that < (7c) < (7s), as one would expect from the defi- 
nition of F (see Eq. (1)) and the fact that (7g) > 0. Now, 
writing (7s) = 1 — (l/(as/3s)), we can estimate the last 
term by assuming (l/(as/3s)) — l/(as/3s). Accordingly, 
we find that on average 



(asps) 



(8) 



((7g) + (7c))' 

1.062 and (7^) > 0.058, 



hence (a^/^s) > 1/(1 - (7g)) 
implying that (7c) > 0-006. 

The linear relation between Pq and Rq is obtained when 
the anomalous years (1991, 1992, 1993, 2001 and 2002, see 
Fig. 1(a)) are excluded from the exponential fit. If these 
anomalous years are included in the fit, an exponential 
regression for Pq (see the thin straight line in Fig. 1(a)) 
yields a non-linear relation between Pq and Rq. In the 
light of these results we may argue that companies show 
typical profit-revenue scenario when both variables are lin- 
early related to each other, at least in an average sense. 
Departures from linearity, as the down-triangles shown 
in Fig. 1(b), may be referred to as extremal, non- typical 
events. Indeed, the cause for such strong deviations from 
the linear relation between Pq and Rq can be traced back 
to specific historical facts^. 

The simple exponential fit for Rq shown in Fig. 1 (a) will 
be used in the following to scale annual revenues and prof- 
its, to take into account the year-to-year variations due to 
the exponential growth in economic activity. One exam- 
ple of such scaled quantities is the revenue itself for which 
we have calculated the probability distribution function 
(PDF) N{r) of scaled revenues r = R/Rq. The results are 
shown in Fig. 2. The scaled revenue PDF N{r) displays 
an intermediate power-law regime with decaying exponent 
~ —1.55 followed by an exponential tail for large r. In 
what follows, scaled profits will be denoted as p = P/Rq. 



^We refer to the gulf war crisis in 1990-1991 and the resulting 
economy's recession, and the post-internct-bubble- and 9/11-effccts 
during 2001-2002. 




Fig. 2: (color online) Distribution function of revenues N{r) vs 
scaled revenue r = R/Ro (circles). The continuous line is a fit 
with the form: n(r) = 0.3 r~'^-^^ exp{~r /6) for r > 0.2. The 
straight line is the power-law form: n{r) ^ (valid within 

the interval 0.3 < r < 2). 



Next, we study the issue of profit fluctuations by 
considering again a generic company having scaled profit 
p and revenue r, at any given year within our database. 
Profit fiuctuations will be considered to be a function of 
scaled revenue r, rather than a function of time. Profit 
fiuctuation, denoted as Ap, is defined as the difference be- 
tween actual p and its expected mean value, here denoted 
as p, i. e. 

Ap = p — p. (9) 

In order to determine the expected mean profit p, we have 
plotted all available values of p and r in our database 
(including the anomalous years) and performed a linear 
regression to the data which should represent the behavior 
of p vs r. The fit, p ~ a + br (not shown here), yields 
a rather small intercept value, a ~ —0.004, which can 
be neglected for our present purposes, while b ~ 0.056. 
Therefore, we postulate that the expected mean profit p 
is a function of r and obeying 



p = (7g)r, with (7g) = 0.052. 



(10) 



We will explain below the reason for choosing the above 
value for (7g). Thus, in the present context, company 
mean profit is a function of solely actual company revenue 
r, times a global market parameter (7g), which is taken 
the same for all companies. Later, we will relax the latter 
condition and discuss the consequences of taken instead a 
different proportionality factor 7g for each individual com- 
pany. Note that by averaging Eq. (10) over all companies 
and years we get (p) = Pq/Rq = (ig), consistent with the 
empirical result Eq. (5). We have also checked that r ~ 1 
and fluctuations of r are essentially stationary all over the 
period considered. 

The fluctuations Eq. (9) can be scaled by using a char- 
acteristic value, such as the standard deviation, Cp, pro- 
vided that the second moment of the distribution be finite. 
Alternatively, one can use a lower-order moment such as 
(|Ap|) to characterize the amplitude of profit fiuctuations. 
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Fig. 3: (color online) (a) Scaled profit fluctuations Ap/^|Ap|) 
vs scaled revenue r = R/Ro, with (lApl) = 0.034. (b) 
Renormalized scaled profit fluctuations A7r/^|A7r|) vs r, where 
Att = Ap/r'', with 77 = 0.6 and (|A7r|) ~ 0.030. 



Values of Ap are plotted in Fig. 3(a) versus scaled rev- 
enue r. As one can see from the plot, profit fluctuations 
are not 'stationary' as a function of r, their amplitudes 
tend to grow with r; the larger the revenue the larger the 
amplitude of profit fluctuations, fn the following, we wish 
to quantify the observed rate of growth of amplitude fluc- 
tuations with revenue and, as a result, being able to find 
out a source of profit fluctuations which is stationary as a 
function of r. To achieve this, we suggest that a suitable 
variable describing fluctuations is given by 



Att 



Ap 



(11) 



with r] > 0. The question arises of how to determine 77. 
To do this, we look at the mean-square fluctuations of 
the data, cr,^ — ((A7r(77))^) - (A7r(7y))^, for fixed 77, and 
search for the minimum of cr^ as a function of rj. Wc find 
a minimum value cr,, ~ 0.05f for rj ~ 0.6. The resulting 
scaled fluctuations Att/ (| Att]) are reported in Fig. 3(b) 
versus r, displaying a satisfactory stationarity. The value 
of 77 thus obtained does not guarantee the vanishing of 
the first moment (Att^. A fine tuning of the value of 
(7g), which enters Eq. (9) and Eq. (fO), can be accurately 
performed in order that (Att^ = 0. This is achieved for 
^7g) ~ 0.052, the value used in the discussions above. 

It is remarkable that 77 7^ 0.5, the latter would indicate 
a standard behavior of fluctuations. The fact that 77 > f/2 
tells us that fluctuations are stronger than one would ex- 
pect if Att were normally distributed (see e.g. [6]). To 
find out the actual shape of the probability distribution 
function, G{e), for e = A7r/(|A7r|) (in the case 77 = 0.6), 
we have plotted it in Fig. 4. As one can see, the shape of 
G(e) is consistent with a power-law decay at the tails with 
a Levy-like exponent a = 1.7. The negative tail of G(e) 
somehow reflects the fact that companies with poor rev- 
enue behavior can be taken out of the Fortune 500 set and 
it may thus become underweighted. Similar arguments 
are used in discussing fund performance (see e.g. [19,20]) 
but a thorough understanding of this phenomenon calls 




Fig. 4: (color online) PDF of scaled profit fiuctuations G(e) 
versus e = A7r/^|A7r|), for 77 = 0.6 (circles). The continuous 
line is a fit with the form; F = a/(l -I- \e/b\ + |e/cp ''), with 
a = 0.55, b = 0.6 and c = 0.9 for e < 0, and F = a/{l + \e/c\'^-'^) 
with a = 0.55 and c = 0.75, for e > 0, implying power-law tails 
with a Levy-like exponent a = 1.7. The dashed line is the 
normal distribution. 



for further studies. To be noted is that processes resulting 
from other human-based activity, such as price variations 
in financial markets (see e.g. [21]), speed fluctuations of an 
ensemble of cars in a closed circuit traffic [22], just to name 
a few examples, also display strongly fluctuating features 
typically characterized by fat-tail distributions. 

We can make contact with the above obtained value of 
77 by arguing that indeed 77 = l/a. To see this, let us write 
the relation e = A7r/(|A7r|) as follows 



P 



(7g)r + (|A7r|) 



(12) 



suggesting that profit p is driven 'deterministically' by the 
first term proportional to revenue r, plus a second fluctuat- 
ing term where amplitude fluctuations are also determined 
by revenue, but to some power 77. From a physical point 
of view, the first term represents a driving or bias field 
and the second one a stochastic part due to an external 
random force acting on the system (see e.g. [23]). Such 
a model resembles very closely the simple approach due 
to Bachelier [24] for describing the temporal evolution of 
stock prices. In our approach, time is replaced by revenue 
and e is not normally distributed (see also [25]). 

Now, imagine we can write the factor r'^e as a sum over 
independent, identically distributed (according to G(e)) 
Levy-like variables e^, such that 



r 'e 



i=l 



(13) 



where a > is a constant and Ur depends on r. A similar 
picture has been used for the description of seed produc- 
tion of forests [26]. Invoking the stability of Levy distri- 
butions, the above sum is also Levy distributed with the 
same exponent a as the single variable, obeying the scal- 
ing relation, J2^=i — ^'nl^" (see e.g. [27]). Even if the 
random variables e are not independent, but long-range 



autocorrelated, the sum scales as 



where the Hurst 
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Fig. 5: (color online) Fraction of profitable (Fortune 500) com- 
panies in the U.S. in the period (1954-2007) (circles) versus 
scaled revenue r = R/Rq. The continuous line is the predic- 
tion from the present work. The dashed line the one from a 
normal distribution of profit fiuctuations. 

exponent H is expected to he H ~ 1/a (see e.g. [28]). 
Details of the corresponding correlations analysis will be 
discussed elsewhere. 

Assuming now that the number is proportional to 
revenue in the form = r/ri,, with arbitrarily small, 
wc find according to Eq. (13), that 

„l/a 

r'^e^ae'^. (14) 

Identifying the parameter a in the form a — r^" and 
noting that e' e, we arrive at the relation i] = 1/a 
claimed above. The values of r] and a obtained here are 
consistent with this prediction. 

In what follows, we will elaborate our findings further 
to consider the issue of profitability, admittedly important 
for being able to make predictions about the probability 
for a generic company to be profitable. This is related to 
the concept of break-even or point at which profits vanish. 
The idea here is to estimate the probability of profitabil- 
ity as a function of revenue r by appropriately taking into 
account profit fluctuations due to statistical market vari- 
ations. 

Our derivation starts from Eq. (12), with i] substituted 
by 1/a, i.e. p = {jg)r + (|A7r|)ri/"e, with a ~ 1.7. At 
break-even (BE), p — and the above relation suggests 
that, e = -(7g)ri-i/"/(|A7r|) = -ebe, where 

which is positive since here (7g) > 0. Now, we define the 
probability for a generic company to be profitable, PpF{r), 
as the fraction of all events for which the fluctuating vari- 
able e > —cbe, since in these cases profit becomes positive, 
p > 0. The probability of profitability can be conveniently 
expressed as the integral (see also [9]) 

/oc 
de G(e), (16) 
"€bE 
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Fig. 6: (color online) Probability of profitability Ppp versus 
revenue [Billions] predicted for 2008 for a typical company. The 
continuous line is the theoretical prediction Eq. (16) and the 
dashed line the linear approximation valid for small revenues 
(here i? < 0.5 B). The dot-dashed line is the prediction from 
the normal distribution and the thin line the corresponding 
asymptotic linear form. 

which is a function of r through esE- One may argue that 
values e — > cxD are not realistic, requiring the introduction 
of an effective upper cut-off for e, that we can denote as 
Ecut- We find that the results become indistinguishable 
from those obtained from Eq. (16) when Ccut > 8 — 10, 
consistent with the range of variations of e obtained in 
Fig. 4. Thus, for simplicity, we take the upper integration 
limit in Eq. (16) as e — » oo. 

Numerical results are shown in Fig. 5, where a compari- 
son is made with the total fraction of profitable companies 
in the period 1954-2007 with respect to the total number 
of companies as a function of scaled revenue. The large 
circles refer to 'typical' years, excluding the anomalous 
ones. If the latter are also included, the fraction of prof- 
itability decreases a bit. Our theoretical prediction works 
satisfactorily well, justifying a posteriori the few ad-hoc 
assumptions made in this work. The dashed line is the 
theoretical prediction in the case in which 77 = 1/2 and 
G(e) be the normal distribution, clearly yielding a poorer 
description of the empirical results. 

Based on these results, we can make predictions for 2008 
using the expected mean revenue i?o — 24.3 B. The results 
of profitability as a function of absolute revenue K (in 
billions) are displayed in Fig. 6. Linear approximations, 
valid for small revenues, are also reported in the plot to 
help making simple estimates for small revenue companies. 

Our results are based on the use of the mean (and 
positive) growth factor (7g) determining expected profits. 
This was done as an attempt to describe a generic, typ- 
ical company. Now, actual companies may behave quite 
differently than this typical behavior. This can be taken 
into account by considering, instead of (7g), the actual 
company-dependent driving factor 7g = 7s — 7c, that is 
the counterpart of Eq. (7) , such that 

P = 7g^ = (Ts-Tc)?^- (17) 
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Here, the factor 7c relates (fixed) costs F to revenue as 
F = 7cii, which in the scaled form becomes f = F/Rq = 

The present break-even and profitability results, valid 
for a typical company, can still be applied to a single com- 
pany with the condition that the break-even value is cal- 
culated according to the particular value of 7g. The result 
is 

(7s -7c)^l-l/a 



(IAttI) 



(18) 



which can become negative if 7s < 7c- In the case 7g > 0, 
also ebe > 0, and the previous conclusions for Ppf(?') 
still apply. When 7g < 0, then cbe < and the lower 
integration limit in Eq. (16) becomes positive, yielding 
lower values of Ppf(?') as compared to Ppf(?') with the 
same revenue r but positive 7g. Thus, for a single company 
the problem reduces to estimate accurately the growth 
factors 7s and 7c in each particular scenario. 

In summary, we have analyzed annual profits and rev- 
enues of U.S. companies over a period of 54 years. We find 
a linear relation between annual mean profit (Pq) and rev- 
enue (Ro), which is at the basis of the concept of typical 
company or mean-profit-revenue relation. In the 'typical 
economy' picture discussed here, expected mean profits 
behave directly proportional to revenue, the latter being 
the driving variable. Furthermore, conjectures allow us to 
study profit fluctuations depending on actual revenue plus 
a fluctuating term governed by a distribution function of 
a Levy type. Strong deviations from the expected typi- 
cal behavior of companies can be referred to as extremal, 
non-typical events. Indeed, within the 54 years data con- 
sidered few extreme events (years) are observed for which 
profits and revenues display strong deviations from the 
linear relation between Pq and Rq. However, such devi- 
ations can be traced back to specific historical facts and 
the corresponding years may be considered as non-typical 
ones. 

Although our conclusions are based on the study of the 

highest revenue companies in the U.S., we believe that our 
results still possess a robust degree of universality to be 
of more general validity, and can be used as a benchmark 
from which one can predict forthcoming profit-revenue sce- 
narios. As an application, the present analysis has been 
used to estimate the probability that a single company 
has in order to become profitable. The analyzed empirical 
data are in support of our suggestions. 



HER would like to thank Claudio Brighenti for illumi- 
nating discussions. 
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